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Abstract 

We introduce the quantitative measures characterizing the rates of decoherence and 
thermahzation of quantum systems. We study the time evolution of these measures in 
the case of a quantum harmonic oscillator whose relaxation is described in the framework 
of the standard master equation, for various initial states (coherent, 'cat', squeezed and 
number). We establish the conditions under which the true decoherence measure can be 
approximated by the linear entropy 1 — Tr/5^. We show that at low temperatures and 
for highly excited initial states the decoherence process consists of three distinct stages 
with quite different time scales. In particular, the 'cat' states preserve 50% of the initial 
coherence for a long time interval which increases logarithmically with increase of the 
initial energy. 
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1 Introduction 



Recently, a significant interest to the decoherence processes in quantum mechanics is observed, 
in particular, due to the problem of stability of quantum superpositions (frequently modeled 
by some kinds of the 'Schrodinger cats' under the influence of the environment p|-p!5 



However, despite that the qualitative picture of the phenomenon seems more or less clear, there 
are no unique quantitative measures of the rate of decoherence or the rate of thermalization. The 
decoherence implies a degradation of the quantum interference effects (manifesting themselves 
in the existence of quantum superpositions) due to the interactions with the 'external world'. 
Since these effects are inherent to the pure quantum states, while they disappear in quantum 
mixtures, it seems natural, on the face of it, to identify the 'degree of decoherence' with the 
degree of 'impurity' of the quantum state, expressed in terms of the 'canonical entropy' [|T^ 
S = — Tr (plnp) or in terms of the 'linear entropy' 0, [l^ s = 1 — Trp^, which is more 



simple for calculations. However, a deeper analysis shows that such an identification leads to 
certain difficulties, especially in the low temperature case. 

Indeed, let us consider the evolution of an initial pure (s(0) = S'(O) = 0) quantum state 
due to a weak interaction with a large reservoir at low temperature. For t > 0, s{t) and S{t) 
assume positive values, so the rate of increase of s(t) or S{t) at t — > can provide us some 
hints to the time scale of the initial phase of the decoherence process [p!T|, H. However, tracing 



the evolution of the entropies for the long time interval, we discover that for a small enough 
temperature of the environment, the entropies, after reaching some maxima, finally decrease to 
very small values which tend to zero when T — > (because the stationary thermal mixed state 
is very close in this case to the ground pure state). Then, identifying the measure of quantum 
impurity with the measure of decoherence, one should accept a strange result that the degree 
of decoherence of the final equilibrium state is almost the same (close to zero) as it was initially, 
despite that the thermal states are usually believed to be the most 'incoherent' 

This example shows that at low temperatures the entropies can serve as the measures of 
decoherence only at the initial stage of the decoherence process. Thus, several questions arise. 
The first one: is it possible to find some other measures which could be used in the whole interval 
< t < oo and for any temperature of the environment? Another question is: under which 
conditions (at which time scale) the usage of the entropies as the measures of the decoherence 
can be justified? 

In the present paper, we answer both questions, introducing the new parameter C whose 
connection with the degree of coherence is indubitable (see section This parameter equals 
one for pure quantum states and for the thermodynamical equilibrium states, for any tem- 
perature T > 0. In section ^ we introduce another parameter V, which can be considered as 
the 'measure of thermalization', since it equals zero for any pure state and 1 for the thermody- 
namically equilibrium state of any quantum system with an equidistant energy spectrum, for 
any temperature T > 0. 

Following the time evolution of the parameters C and V in the process of the thermal 
relaxation of various initial states (Fock's, coherent, squeezed, 'cat') of the harmonic oscillator, 
described in the framework of the 'standard master equation' (sections ^ and ^), we find the 
conditions, under which the linear entropy can serve as a reasonable 'measure of decoherence'. 
Moreover, in section ^ we demonstrate that the decoherence of highly excited initial states at low 
temperatures goes through three distinct stages, characterized not by some unique 'decoherence 
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time', but at least by two times with quite different dependences on the initial energy and the 
temperature. 

The first time ti (which is usually identified with the time of decoherence) is, roughly 
speaking, inversly proportional to the product of the energy of quantum fluctuations by the 
number of photons per mode of the reservoir. During this short time interval the parameter 
C rapidly decreases from 1 to some finite value which depends on the initial state. Then 
C{t) remains at a more or less constant level (for the 'cat' states) or even can increase with 
time (for the 'squeezed' states). And only after the 'ultimate decoherence time' ta, which 
increases (logarithmically) with the increase of the initial energy, the coherence coefficient goes 
monotonously to the final zero value. 

The evolution of the 'thermalization parameter' I'(t) is similar to certain extent to the 
behaviour of 1 — C{t): the initial rapid increase from zero to some intermediate value, then 
some stabilization or even decrease, and the final rapid transition to the equilibrium unit value 
after the 'thermalization time' tx, which also depends logarithmically on the initial energy. The 
difference between td and consists in their temperature dependences: t^ has a finite limit 
value when the temperature T tends to zero, whereas is inversly proportional to the absolute 
temperature at T — >■ 0, in accordance with the third law of thermodynamics (the inattainability 
of the absolute zero of temperature implies that the rate of the relaxation processes must go to 
zero as T ^ 0). 

2 The measure of 'coherence' 

The controversies of the identification the decoherence measure with the von Neumann' or 
linear entropies take their origin in the invariance of these entropies with respect to the choice 
of the basis in the Hilbert space: the entropies do not distinguish the equilibrium state (or 
other stationary states in the case of 'colored' or 'squeezed' reservoirs) from any other mixed 
one. But such a symmetry with respect to the choice of the basis in the Hilbert space is 
obviously broken in the relaxation processses, when all possible initial states tend to the unique 
equilibrium state, whose density matrix is diagonal in the distinct basis formed by the energy 
eigenstate of the Hamiltonian (or by some other distinct basis in the case of more sophisticated 
artificial reservoirs considered recently). Thus it seems natural to suppose that the measure 
of decoherence must depend explicitly on this distinct basis (the concept of the broken 

symmetry of the Hilbert space was used as a basis for introducing the polarized distance between 
different quantum states in [0). In some special cases, when it is known exactly to which 
specific family of quantum states (e.g., coherent state, even/odd coherent state, etc.) the 
initial quantum state belongs, the prefered basis may be different from the energy one, so 
that some special approaches can be used, as well. We shall discuss such a situation later on. 
However, in the generic case, when the type of the quantum state is not known beforehand, the 
only available information is contained in the set of matrix elements of the statistical operator 
with respect to the physically distinguished energy basis: p{t) = J2mnPmn\fn){n\- Therefore, 
it seems reasonable to define the measure of coherence in terms of the coefficients pmn- Since 
the decoherence is usually identified with the disappearence of the off-diagonal elements of the 
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density matrix ||pmn||, it is natural to define the normalized measure of coherence as 

C{t)= E Kni'W/E \pmnf{0). (1) 

Then C(0) = 1, while C = for any 'completely incoherent' state without off-diagonal matrix 
elements in the energy basis (provided initially at least one off-diagonal element was different 
from zero). Introducing the 'diagonal part' of the operator p 

= EP"I^)("'I' Pn={n\p\n), (2) 

n 

and taking into account the property Tr [ppd) = Tr (p^) one can rewrite ([^) in the form 

^ ^ TT[p{t)-pd{t)? ^ pit) - A(t) 
Tr[p(0)-p,(0)]' MO)-A(O) 

p^Trp^ \^i:vpl = Y.A (4) 

We shall call p the 'total purity' and A the 'diagonal purity'. In many cases of practical interest 
both 'purities' can be calculated rather easily. For example, if one knows the Wigner function 
ni EDI (we assume h = 1) 



Q+l) (5) 



W{q,p) = J dv e*^^(g-| 

Trp = J W{q,p)dqdp/ {27r) = 1 

p = J W\q,p)dqdp/{2'K). (6) 

As to the 'diagonal purity', it can be calculated either by means of a direct summation of the 
series in (§), or, equivalently, it can be expressed as the integral 



then 



/o 27r 
where 

oo 
n=0 

is the diagonal generating function. Below we use both methods. 



3 The measure of thermalization 

A qualitative measure of thermalization can be introduced in the following way. The analysis 
of the low temperature behaviour of the entropies shows that the troubles mentioned in the 
introduction arise due to the double nature of the ground state, described by the density 
operator po = |0)(0|. On one hand, this state is pure, with Trpg = 1. On the other hand, it is 
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the limit of the equilibrium states, which are conceived to be completely decoherent. Therefore 
it seems reasonable to exclude the state po in some way. One of the possibilities is to take 
a simple expression for the linear entropy and to divide it by a proper time-dependent factor 
which would ensure a nonzero limit at t oo. This goal can be achieved, for instance, if 
one chooses as the normalizing factor the Hilbert-Schmidt distance between the states p{t) and 
Pq. If the system under study has a finite number of energy levels (e.g., spin systems), then 
there are some grounds to treat the state with the maximal energy Pf = \Emax) {Emax\ on the 
same footing as po- Thus we arrive at the parameter (introduced for the first time in ETI, but 
identified erroneously with the measure of decoherence) 

1 - Trp2 
V = 



1/2 

Ti{p- PoYTiip- PfY' 



'-"^ (9) 



M 



where po = Tr (ppo) = (0|/)|0) is the ground state occupation probability, while pj = Tr(/3pj) 
is the occupation probability of the level with the maximal energy (evidently, for quantum 
systems with infinite dimensional Hilbert spaces, such as a harmonic oscillator, = for any 
physical state possessing finite energy). 

For systems with equidistant spectra, En+i — En = AE = const, the equilibrium occupation 
probabilities read 

Pn = ai-e)/(i-r) 

where M is the total number of levels, n = 0, 1, . . . , M — 1 and ^ = exp{—(3AE) < 1 is the 
Boltzmann factor. Then 

/^e,= (l-0 (l + r)/[(l + (l-e 

and we see that V^q = 1 for any value < ^ < 1, while X> = for all pure states. For 
this reason, we may consider V as 'the measure of thermalization'. Note that we have an 
indeterminacy in (P) if ^ = T = 0, thus the case T = must be excluded. But as we know the 
limit of exact zero temperature is an idealization, thus we do not have to worry on this issue. 

For systems with nonequidistant spectra the value of V^g depends on temperature, nonethe- 
less the limits at T — and T — cx) still equal 1. For T — oo we have po = pi = . . . = Pf = 
1/M, consequently p = M ■ (1/M)^ = 1/M and T> = 1. In the low temperature case T — 0, the 
equilibrium statistical operator is close to po|0)(0| with pi ^ 1 (where |1) is the first 

excited state), while the contribution of other states can be neglected (we consider the systems 
with discrete energy spectra). Then po+pi = I, p = p^+pf, 1 + p — 2po = 2pl, 1 + p — 2pf = 2, 
and 1 — p = 2pi (up to higher order terms). As a result, we have V = 1 again. 
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4 Decoherence dynamics 



4.1 Time evolution of the quantum state 

We confine ourselves to the analysis of the time dependence of the 'coherence coefficient' C 
in the process of thermal relaxation of the harmonic oscillator described in the framework of 
the standard master equation ||2l], ^ (more general models were considered, e.g., in 



(10) 



dp/dt = 7 (1 + z/) (2dpd^ — d^dp — pd^d^ 
+7Z/ (^2d^pd — 



dd^ p 



paa^ 



a'a, p 



Here d and d'^ are the usual bosonic annihilation and creation operators, v is the equilibrium 
mean number of quanta in the reservoir corresponding to the given mode, and 7 > is a 
damping coefficient {h = u = 1). 

An immediate consequence of equation (p!OD is the universal expression for the purity loss 
rate in an initial pure state p^ = p = \'ip){'ip\ (cf. JTI], |l^ ^ ^] ) 



p 



2Tr (pp) ^^^ = -47 [z/+(l + 2z/)orJ, 



'111 



where cxa = (^Vl^^^lV') ~ I (V^I^IV') 1^- The 'primary' purity loss rate is minimal for the coherent 
states with aa = 0. In the generic case it is roughly proportional to the average number of 
thermal photons in the reservoir and to the 'reduced' energy of quantum fluctuations in the 
initial pure state, 



£0 



cr, 



(0) 



a 



(0) 



(12) 



where cr^°^ and a^^ are the variances of the quadrature components q = (a + ) /V2 and 

p = i (^d^ — a) /V2 in the initial pure state 

Due to equation (p!TD the initial evolution of the 'purity' has the linear dependence on time 
p{t) = 1 — t/ti + ■ ■ ■, where the 'primary purity loss time' equals 

t, = (47)-i[z/ + (l + 2z/)aJ-\ (13) 

Note that some 'microscopic' models, based on an explicit coupling of the system under study 
with a large reservoir, result in a quadratic time dependence p{t) at t ^ fl^. This apparent 
contradiction is explained by the fact that the 'microscopic' models and the phenomenological 
master equations describe the evolution of the subsystem in different time scales. Actually, the 
master equation describes a 'coarse-grained' evolution averaged over many periods of the fast 
oscillation, so the physical meaning of the limit t ^ in the case of the master equation is 
quite different from the same (formally) limit in the 'microscopic' models. 

To calculate the time dependence of the 'purity' p in the whole interval < t < cxo with 
the aid of formula (|^) we need the time dependent Wigner function W{q,p,t). It obeys the 
Fokker-Planck equation which follows immediately from (|10[): 



dW 



— ([7g - p\W) + Q^{bP + iW) 

Qq2 Qp2 j ■ 



+7 + 



(14) 
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The solution to equation (0) can be written as 

W{q,p;t) = J IC{q,p;t\q',p',0)W{q',p';0)dq'dp'. 



(15) 



The propagator }C{q, p; t\q', p', 0) was calculated by means of different methods in ^ ; 
the explicit expressions in the case of most general multidimensional time-dependent quadratic 
operator (with respect to q, p, d/dq, 8/ dp) in the right-hand side of the Fokker-Planck equation 
were given in In the case involved the general form found in is reduced to (see 

appendix) 



/C(g,p; 0) = (27rcr*M) ^ exp | — (2cr*M) ^ 



2 I 2 
% +Pt 



(16) 



where a^ = v+\, 



qt = qcost — psint, pt = pcost + qsint, 
and the 'compact time' u is given by 

u{t) = 1 - e-^^K 



(17) 



The consequence of the master equation (|TOD is the closed set of equations for the diagonal 
elements of the density matrix in the Fock (energy) basis 

Pn = 27(1 + U) [{n + l)Pn+l - npn] 

+2-fu[npn^i-{n + l)pn]. (18) 

These equations, in turn, are equivalent to the simple first order partial differential equation 
for the diagonal generating function (P) 



^ = 27(1 -z)[l + v{l - z)]^ - 27z/(l - z)G 



(19) 



The solution to (|19|) reads |^T], ^ ^ 



-Gn 



' z + u{l + - z) 



(20) 



1 + uu{l - z)~" \ 1 + uu{l - z) 
where Gq{z) = G{z, 0). Putting 2; = 1 in ( pO|) we verify the normalization condition G(l, t) = 1. 



4.2 Initial coherent states 



As the first example we consider the evolution of the initial coherent state \a), a 
Applying the propagator ( PIB] ) to the initial Wigner function 



aexp[t(f 



iy('=°'^)(0) = 2 exp \-{q- V2^cos(f)f - (p- v^sin, 
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we obtain 



W^'°^\q, p, t) = exp - g - v^e"^* cos(0 - t) 



p — v2ae ^* sin(0 — t) 



where 



(21) 



(22) 



Uu) = {l + 2uuy^ 

coincides with the 'total purity': = ^„{u). 

To calculate the 'diagonal purity' we use the explicit expression for the time- dependent 
diagonal matrix elements in terms of the Laguerre polynomials pl| , |33[] 



p 



(coh) 



\n+l 



exp 



a{u — 1) 
1 + uu 



a{u — 1) 
uvil + up) 



(23) 



Then the sum in (^) is reduced to the known series 134 

oo 

Y^Ln{x)Ln{y)z- 



l-z) 



exp 



n=0 



X + y 



'l-z 



(24) 



(Jo(z) is the modified Bessel function), so we obtain 

A(™'')=e.(w)exp(-r7)/o(r/), (25) 

r] = 2a{l-u)Uu). (26) 

For a ^ 1 it is sufficient to take into account the first terms of the Taylor expansion of the 
function e~^Io{ri) to obtain 

C(^°'')(t)^e^H(l-^)=^'We-'^*. (27) 

In this case the time dependence of the 'purity' has a little in common with the time dependence 
of the 'coherence'; the same is true even for a ~ 1: see figure ^ 

For the highly excited initial states with a ^ 1 the asymptotics of the modified Bessel 
function, /o(a;) ~ {2Trx)~^^'^e^ for a; 3> 1, yields 



1 



2^/a7^ 



VI + 2Mz/e^* - 1 



In this case, the contribution of the diagonal elements to the total purity is small compared with 
the contribution of the off-diagonal terms, therefore the sum over m 7^ n in (|ip can be replaced 
by the sum over all values of m and n, and the correlation coefficient can be approximated by 
the purity /i. This is just the case considered in most of the papers devoted to the decoherence 
of initial 'macroscopic' quantum states. However, C ~ /i only under the condition A <^ /i. 
Consequently, the linear entropy 1 — /i can be considered as a measure of decoherence only for 
a 3> 1, and under the additional restriction a^,^(M) exp(— 27^) ^ 1, i.e., at the time scale 



t < ^ (27)-! ln[(a + 2z/)/(l + 2z/)]. 



(28) 
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For t > tt, the identification of tlie 'purity' with the 'coherence' leads to incorrect results as 
mentioned in the introduction. 

The 'final decoherence time' td can be defined by means of the equation 

C{td) = Pfieg (29) 

where < 1 is some small number whose choice is a matter of convention (say, = 0.1), 
and fieq = (1 + 2z/)~^ is the equilibrium value of the 'purity' (we take into account that C is 
proportional to = fi^^"^\ according to equation (pSD ). To solve equation ( ^9]) for sufficiently 
small f3 we can use the asymptotical form of the coherence function at 7^ ^ 1 (when u ^ 1) 

&°^\t) ^ 2a(l + 2v)-^e-^^\ a > 1, 7^ > 1. 

Thus we obtain the estimation 

^ (27)-Mn[2a/i,,//?] (30) 

which holds for a^^q 3> 1. The evolution of the 'total purity' and the 'coherence coefficient' for 
highly excited initial coherent states at different temperatures is shown in figure 0. 



4.3 Initial 'cat' states 

Now let us consider the family of the initial 'Schrodinger cat' states [H, ^ § 



■^)=Ui\a) + 



a) 



a = \a\ 



(31) 
(32) 



AT = (2 [1 + cosy?exp(-2a)])~^/^ 

The special cases of this family are even (yj = 0) and odd ((/? = tt) coherent states and the 
Yurke-Stoler states (v? = 7r/2) 0. The Wigner function of the state (|3T|) reads (hereafter we 
assume a = -y/a to be real) 



W^^"*^ (g, p; 0) = 4Af^ exp (-g^ _ 



X 



e " cosh(v8ag) + cos(v 8ap + </)) 



(33) 



Applying the propagator (|16|) to this function we obtain (see also |]T0| 

iy(^"*)(g,p; t) = 4Af%exp [-^, (<f + 
X |exp [— 2a(l — u)^y\ cosh 

+ exp [— 2a(l + 2u)u^y\ cos 



iyJSail - u) qt 



-u)pt + ip 



(34) 



where the function C,u{u) was defined in equation (|2^) . Calculating the integral in (P) we find 
the 'total purity' 



= 2^^%{u) [1 + 4cos(^e-2'' + cos(2^)e"^" 
+ exp [-4a(l - u)C,u] + exp [-4aM(l + 2z/)^^] 



(35) 
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The photon distribution function can be written as 

p^*) = [p^:''^\a) + cos^e-'y^:^^\-a)\ , (36) 

where p^n°^\a) is given by (P^. Calculating again the sum I]Pn "with the aid of formula ( |2^ ) 
we find 

y{cat) ^ 4Ar4^^(M)|/o(r/) [e-^ + cosVe''"^'' 

+2cos(/?e"2Vo(r/)| (37) 

where Jo (-2) means the usual Bessel function, and r] was defined in equation (pG]). If a <^ 1, 
then we have ~ ^!y(^^) and the same formula (pTf ) for A. 

For highly excited ('macroscopic': a ^ 1) initial cat states the phase ip becomes unimpor- 
tant, and A/"^ ^ 1/2. Until a(l — -u) = aexp(— 27t) ^ 1, the 'diagonal purity' is small, similarly 
to the case of coherent states: A'-^"*^ ~ exp(7t)/Y^. Then the coherence coefficient can be 
replaced by the 'total purity' /i, which rapidly decreases from 1 to the value \iuiu): 



~ + exp [-Aau{l + 2u)i^{u)] 



staying at this level until a{l — u) becomes smaller than 1. In particular, in the low temperature 
case 1/ -C 1 we observe the 'plateau' /i ~ C ~ |: see figure (|^). 

At the final stage of the evolution n goes to the equilibrium value = (1 + 2z/)~^ as 

~ \iu{u) (1 + exp[-4a(l - u)^y{u)] 



2 

but now it becomes compatible with the 'diagonal purity' A ~ iu{u) exp[— 2a(l — u)^y{^i)], so 



C^^-') ^ \Uu){l - exp[-2a(l - u)Un)]Y ■ 



For a(l — m) ^ 1 we have C*^™*^ ~ 2a^/igg exp(— 47^). Then equation (|29|) yields the 'ultimate 
decoherence time' t^'^"*^ ~ (27)~Mn a^eqJ^j^ (if a/igg ^ 1), which is only slightly less than 



the similar time t^^°^^ (|3 



4.4 Decoherence in the accompanying basis 



In the special case of the 'cat' states there exists another natural choice of the 'diagonal' part of 
the statistical operator, different from (0). Namely, one can define the 'accompanying' diagonal 
operator as pmix = | (|ct)(«| + I ~ a){—a\). The Wigner function of this quantum mixture is 
the sum of two coherent Wigner functions (0) with opposite values of parameter a: 

1 



= 2^^{u) exp {~^^{u) [g2 + p2 + 2a(l - u)] } 



X cosh 



^,y{u)^j8a{l - u) qt 



(3^ 
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Then the 'accompanying' (normahzed) coherence coefficient can be defined as 



X 



(39) 



Performing the calculations we obtain 



i^[u) (1 - exp [-4a(l - u)iJ\) 



cos^ V9e ^'*) 



(exp [— 4au(l + — cos^ y^e' 



-4a 



(40) 



The time evolution of this parameter is essentially different from the evolution of the coefficient 
C: see figure (^). This is not surprising, since the phenomena observed from the moving 
('accompanying') frame of reference in many cases look very different, compared with their 
appearance in the fixed frame. 



4.5 Initial squeezed states 

Now let us consider the states possessing the Gaussian Wigner functions 29, 35, 3C 



W{q, p) = d exp I - ^ [ap(g - qf 



-2aqp{q -q){p-p) + (Tgip - pf 



(41) 



the (co)variances of the quadrature components, whereas q and 

al must satisfy 



where a^, and cr^p — ^pq 

p stand for the average values of these quadratures. The parameter d = (jpG^ 



, ^ d>l/A. It is related to the 'purity' 



the Schrodinger-Robertson uncertainty relation 
of the state as /i = (4(i)~-^/^. In the thermal state, a^"^^) = a^'''> = = ^ + u, apg = q = p = 0. 

The evolution of the five parameters determining the Gaussian state is governed by the set 
of equations following from the master equation ( [T(]| ) 



dq/ dt = p — 7g, dp/ dt = —q — •yp 



(42) 



The solutions read |p9| 



cr, 



qp 



m 
Pit) 



2(T,p-27(T, + 7(l + 2z/) 
2cT„„ - 27ap + 7(1 + 2z/) 
- 270-5^ 



cr. 



'qp 
On 



e "^^ [go cos t + pq sin t] 
e~^* [po cos t — qo sin t] 



(43) 
(44) 
(45) 

(46) 
(47) 



a,(t) = a, + e-2^*[(aW-a,) 



a^"^ -cr*)sin't + cr(5)sin(2t) 



(4J 
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cTpit) = a, + e-^"" [(af^ - a,) cosH 



= e-^^* [aW cos(2t) + - (a^ - a^) sin(2t) 



(49) 



(50) 



Any initial pure Gaussian state is unitarily equivalent to the squeezed state, defined as an 
eigenstate of the canonically transformed operator b = cosh pa + sinh p with a complex 
eigenvalue a = y/aexp{i(f)) and a real 'squeezing parameter' p (sometimes it is called also a 
'two-photon state' [051). Therefore we parametrize the initial variances and average values as 



< = 0, 



2 ' P 2 

go = V2ae~'' cos 0, po = V2ae^ sin (j). 
The 'total purity' p does not depend on the first order moments of the coordinates: 

'1/2 



p 



(1 + 2w)^ + 4m(1 -m)(1 + 2z/) sinh> 



(51) 



On the contrary, the 'diagonal purity' A depends on all 5 parameters of the one-dimensional 
Gaussian Wigner function. 



The generic diagonal generating function found in can be expressed as 



D 



zgi - z^g2 



90 



where 



;i + zf + Ad{l - zf + 2(org + ap) (l - z^) 



D = 1 + 2(a„ + a„ 



Ad 



go = p\2ag + 1) + q\2(rp + 1) - Apqa. 



(52) 



pq 



91 



2p2 



+ 2f 



+ ^lq + ^P+^ 



-Apqapq{aq + (Xp + 1) 



92 = 2p''{al + alq--^+2q^{al + alq-\ 



-Apqapq{aq + dp) 
For the initial pure squeezed states G{z,t) can be written as 

2\-V2 



G{z; t) = (^f -hz + cz- 

F~Bz + Cz^ 
f — bz + cz"^ 



X exp 



—ail — u) 



(53) 
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where 

/ = + + (1 + M + 2w) sinhV, 

b = 2uu{l + uu) + 2u{l + 2z/)(l - u) sinh^, 
c = (uuy — (1 — -u) (1 — M — 2uu) sinh^ p, 

F =^[l-u + R{l+u + 2uu)] 

B =l-u + Ru{l + 2v) 

C = R{l-u-2uu)] 

R = cosh(2p) - sinh(2p) cos(20) 
Using formula (|^) and the relations 

/-6 + c=l, F-B + C = 

we obtain after some algebra the following integral representation for the time dependent 'di- 
agonal purity' of a generic initially squeezed state: 



A 



2^7 



n^^) 



exp 



-4a(l -m) (V + Ysm'^j 



sin^ 7 



(54) 



where 



$(7) = (l + 2b sin^ 7] + (1 - m)2 sinh^(2p) sin2(27) 
V = R{1 + 2uu) + {1-u)(R-1 + 4i?sinhV) 



Y 



AuB 



z/(l + Mz/) + (1 - m) (1 + 2z/) sinh^ p 



+2(1 -m) (1 -i?-2i?sinh^ 



P 



The evolution of the 'purity' and the 'coherence coefficient' for the highly squeezed (p > 1) 
initial state is illustrated in figure (|]). In this case p ~ C up to the values of the dimensionless 
time r = 27^ ~ 1. In contradistinction to the cases of coherent or 'cat' states, the coherence 
coefficient is not monotonous function of time, but it tries to follow the increase of the 'purity' 
at r > 1, before going finally to zero. 

The integral ( |5^ ) can be easily calculated in the long-time limit l — u = exp(— 27t) = e <^ 1 
at zero temperature (z/ = 0): A = 1 — 2e {smh? p + aR^ + 0{e'^). Comparing this expression 



with the similar expansion of the 'total purity' p 
where 



)1|) we obtain C ^ p — X ^ 2eaR = 2e£a 



aR 



(55) 



is the initial 'classical' energy (the total energy without the contribution of the vacuum fiuc- 
tuations). Then equation (|29|) yields the 'ultimate decoherence time' t^^'^^ ~ {2'y)~^ In {2Sci / (3) 
which has the same order of magnitude as the time t^^"'^^ (^) for the coherent state with the 
same 'classical energy' |ap (at zero temperature peq = 1)- If a = (the initial squeezed vacuum 
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state), then one should calculate fi and A up to the second order terms with respect to e. In 



this case we obtain C ^ fi — X ^ ^e"^ sinh^(2p), and t^^''^ ~ (27) ^ In sinh(2p)/2v^ 



Since we 



consider the case p ^ 1, we can replace sinh(2p)/2 by sinh (p) = E — ^ ^ 
total energy in the case discussed. 



E, where i^^ is the 



5 Thermalization dynamics 



The concrete evolution of the thermalization coefficient depends on a peculiar 'competition' 
between the 'total purity' p(m) (which was calculated in the preceding sections) and the ground 
state probability po{u), which can can be easily found from the 'diagonal generating function' 



Po{t) 



1 + uu{t) 

In the case of initial coherent state we have (a 



Go 



'uit)il + u) 
1 + uu{t) 



lap) a 



(coh) I 



(coh) 

Po 



1 + uu 



exp 



a(l — u) 
1 + uv 



so the 'thermalization coefficient' reads 



l + 2nz/ / 
1 + — — — 1 - exp 



[uv] 



a(l — u) 
1 + uv 



(56) 

exp [a{z — 1)]. Consequently, 

(57) 

^ -1/2 



Evidently, the case i/ = should be excluded in this expression. The dependence on the 
displacement parameter a disappears for a(l — m) ^ 1, when all the functions T)^"^\u) merge 
to 'D^^\u) ~ uv/{l + uv). This slow evolution is transformed into a fast transition to the 
equilibrium value if a(l — m) ^ 1: 



For the 'cat' states equations ([ 



1 + 2v as 
2(1 + z/)^ 

and (|36D yield 



(5^ 



(cat) 

Po 



2M^ 



1 + uv 
X < 1 + cos exp 



exp 



a(l — u) 
1 + uv 
2au{l + v) 
1 + uv 



(59) 



At low temperatures {v <^ 1) all the exponential functions 'die out' if a(l — ■u)/(l + z/) ^ 1 and 
au/{l + z/) ^ 1, and we observe the 'plateau' 'D^'^°-^\u) ~ | (see figure H). 
For the squeezed state, formula ([53| ) yields 



Po 



G{0,t) = r'/'exp [-a{l-u)E/f]. 



(60) 



For large values of the squeezing parameter p and v <^ 1, the 'purity' (which does not depend 
on a) is given by p ~ (2 sinh p\Ju{l — u)^ <^ 1, unless u is close enough to or 1. If 
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a{l — u) 3> 1, then po 1, and we observe a universal (independent of a) behavior of the 
thermahzation coefficient 

pN)/ A 1 -/^ 4M(l-M)sinhV 



1 + /i 



1 + 4m(1 - u) sinhV + 1 



The P-factor rapidly increases for a small time interval t < ti ~ (7sinh p)^^ (note that 
sinh^p is just the 'reduced' energy of fluctuations in the initial squeezed state, Sq — For 
t > ti we observe some 'plateau', whose extension corresponds approximately to the interval 
0.07 < u < 0.93 (inside this interval, the values of the function u{l — u) are not less than a half 
of the maximum value at m = 0.5). For the values of u close to 1, the "D-factor may decrease, 
following the decreasing linear entropy, but finally the term po enters the game and prevents 
the thermahzation coefficient from falling down to zero. This final stage of evolution seems 
very fast in terms of the 'compact time' m, but it is not so dramatic with respect to the scaled 
time r = 27/!:: see figure |^. 

It is interesting to consider also the thermahzation of the initial M-photon Fock state |M). 
In this case the off-diagonal elements of the statistical operator in the Fock basis are equal to 
zero identically for any time t > 0, so the 'total purity' p coincides with the diagonal one A. 
The initial diagonal generating function equals Gq{z) = z^'^ . Consequently, 

[l + uu{t)f'^^ 

whereas the integral (|^) can be transformed to the form 

lo 27T [c — d COS ip] 

a = u^{l + z/)^ + (1 - M - uuf, d = 2uz/(l + uu) 
b = 2m(1 — u — uu), c = v^v? + (1 + uv^'^ 
To calculate the integral ( |6TD we designate it as Im and introduce a new generating function 

dip 



Aw=rg .'''+.''^°^i,. (61) 



Q{y) = y: W = / 



2'n-[c — ya — {d + yb) cos ip] 



The last integral is given by the expression [34 



Q(y)=[{c-yay-{d + yby] 

which has the same structure as the known generating function of the Legendre polynomials, 
so after some algebra we obtain 

\l-2u{l + iy)\'\ ( {l-uf + u\l + 2uf \ 



(1 + 2mz/)^^+i \{l + 2uu)\l-2u{l + u)\J 

The typical dependences 'D(u) for the coherent and Fock states with different initial energies 
are given in figure ^. 
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6 Three stages of decoherence and thermalization 



We see that at low temperatures the decoherence and thermalization of highly excited initial 
states go through three distinct stages. The first one is rather short, its characteristic time being 
determined completely by the initial energy of quantum fluctuations, ti ~ {''yS)~^. However, 
the coefficients C and V do not assume their equilibrium values (0 and 1, respectively) at the 
end of this stage, but they remain approximately constant for a rather long period of time. The 
total destruction of coherence is observed only after the time t^ ~ (27)"^ In (E) ^ ti, where E 
is either the total energy or its 'classical' part (depending on the initial state). This time tends 
to a finite limit when the temperature T goes to zero. 

The disappearance of the off-diagonal matrix elements of the statistical operator (decoher- 
ence) does not mean that the energy level populations reach their equilibrium values. This 
happens only after the 'thermalization time' t^, which can be evaluated from the asymptotical 
behavior of the thermalization parameter P at t ^ 00 in the form of the Taylor expansion with 
respect to the small variable e = 1 — u = exp(— 27/^). For example, for a generic squeezed state 
with nonzero mean values of the quadrature components we have (see also 



where the 'classical energy' is given by (|55|). Assuming (for z/ ^ 1) {Sd/i^"^) exp(— 27^) ~ 1 we 
obtain the estimation tr ~ (27)"^ In (£^c//^^) which shows that the 'thermalization time' may 
exceed essentially not only the decay time 7"^ but the 'ultimate decoherence time' td, too. In 
particular, for the initial coherent state {p = 0) we have tr ~ 7"^ In 

The situation resembles the classical theory of magnetic relaxation, where we have also two 
characteristic times: the time of transverse relaxation (dephasing) T2 (analog of t^) and the 
time of longitudinal relaxation Ti (analog of tx)- The difference is that in our case both times 
depend not only on the properties of the environment (through the constants 7 and z/), but 
also on the initial state (through its energy). Besides, there exists the third time - the 'primary 
decoherence time' ti. 

For the states with zero mean values of the quadratures, the expansion of 1 — P begins with 
the quadratic term = exp(— 47t). For the initial vacuum squeezed state {a = 0) we have 



D = 1 - 



sinh(2p) 



+ O [e^) . (63) 



4z/(l + i^) 

If p > 1 and z/ ^ 1, one can rewrite ( |63D as 

V^l- [Ee/{2v)f (64) 

where E = sinh^ P + I ~ i 6^P(2p) is the total energy of the initial state (it coincides with the 
energy of fiuctuations in the case involved). Consequently, fr^}''""'^ ~ (27)"^ In (E/u). 

A similar behavior of T^lu) at 1 — n ^ 1 is observed for the 'cat' states. If a ^ 1, then the 
dependence on the phase if becomes unimportant, and we obtain 

V = 1- - — — — ■ + 

4z/2(l + Z/)2 
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In this case the total energy E ^ a ^ 1, and we arrive again at the equation (0) (if z/ ^ 1) 
which yields ~ (27)"^ In (E/u), similarly to the case of the vacuum squeezed state. 
For the Fock states we obtain 



If M 1, then E ^ M, and we have t^"'^'' ~ 7"^ ln(E/-y/z/). We see that the 'thermalization 
time' tx depends logarithmically on the initial energy. Besides, it has the strong tempera- 
ture dependence, growing as at T — > (remember that u = [exp{huj/kBT) — so 
u ^ exp{—huj/kBT) at T ^ 0), in a complete agreement with the third law of thermodynam- 
ics. The dependence of the 'thermalization time' on the mean equilibrium photon number u 
enables ordering different families of quantum states with respect to their robustness against 
the thermalization (while the 'primary time' ti is the same for all states with equal values 
of the energy of quantum fluctuations). The coherent states are the most robust ones, then 
follow squeezed and 'cat' states, whereas the Fock states, being 'the most unclassical states', 
are thermalized much faster than all the others. 

Another interesting feature of the decoherence and thermalization process is the existence 
of 'plateaus' in the dependences C{u) and T>{u) for several different types of states (excluding 
the coherent states) possessing high initial energy (provided the temperature is low enough). 
In the cases of the squeezed and Fock states the altitudes of 'plateaus' tend to 1 for V^u) and 
to for C{u) when the initial energy increases. But for the coherent 'cat' states, the metastable 
values of the coherence and thermalization coefficients remain finite even for E ^ 00: Cpu ~ \ 
and Vpit ~ |. Consequently, some degree of coherence (with respect to the fixed energy basis) 
survives in the 'cat' states for a long period of time t < td ^ (27)~Mna. Perhaps, this fact 
could be important for applications. 



7 Conclusion 

We may conclude that the new quantitative measures of decoherence and thermalization shed 
new light on the details of the decoherence process accompanying the 'standard' thermal re- 
laxation of a quantum harmonic oscillator, showing that this process has three distinct stages 
in the case of highly excited initial pure states and low temperatures. In particular, we have 
shown that at low temperatures the 'ultimate decoherence' is achieved after rather long inter- 
val of time, which is essentially greater than the relaxation time and the 'primary decoherence 
time' which was the central subject of previous studies. Our analysis permits to find the con- 
ditions under which the 'purity' or the 'linear entropy' can serve as reasonable measures of 
(de) coherence: the initial energy of the quantum state must be much greater then the mean 
energy of the reservoir oscillators, Eq ^ l + 2z/. However, even under this condition the 'purity' 
can be used only to describe the initial stage of the relaxation process, but it cannot replace 
the true measures of 'coherence' for the whole time interval. 
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A Propagator of the Fokker-Planck equation 



Under certain conditions the process of relaxation of linear r- dimensional quantum systems 
(such as a system of coupled oscillators or a charged particle in a homogeneous electromagnetic 
field and in a confining parabolic potential) can be described in the framework of the Fokker- 
Planck equation for the Wigner function ||2^, P2| -P5| 

dW d d'^W 

^\W]+D,,^^ (A.l) 



dt % ' dvidyj 

where = l,2,...,2r; the 2r-dimensional vector y consists of the linear combinations of 
the Cartesian coordinates qi and canonically conjugated momenta Pi (in the simplest case 
y = (ct, p))- The drift matrix A and vector K do not depend on the phase space vector variable 
y, although they may have, in general, arbitrary dependences on time. However, the diffusion 
symmetric matrix D = ||-Di,|| cannot be arbitrary, since the physically acceptable solutions to 



equation ( [A.l| ) must satisfy the condition of the positive semidefiniteness of the corresponding 
statistical operator. This condition is fulfilled provided the matrix D* = D + ^ ^AS + EA^ 
is positively semidefinite ^ The elements of the antisymmetric c-number matrix 

E = ||Ejfc|| are the commutators Sj^ = f [yj,yk\- In the case of the single space coordinate the 



matrix condition D^, > is equivalent to three scalar conditions ^6] , ^ |^ 

16 
D 



DppDqq 



D^^^detD> 



(A.2) 



D„ > 0, 



> 0, 



D 



D 



pq 



^pq ^qq 



Since ( |A.1|) can be considered as the Schrodinger equation with an effective quadratic (although 
non-Hermitian) Hamiltonian, the propagator G'(y,y',t), 



W{y,t)= Jdy'Giy,y',t)w(y',0), 



can be found with the aid of the method of quantum time- dependent invariants given in [ pO| , 
However, to find its explicit expression it is sufficient to know that this propagator is 
a Gaussian, so, as any Gaussian Wigner function [EOl BO, |5^, |5^ it can be written as 



(A.3) 



G'(y,y',t) = (27r)-'^ [detMM''^' 

X exp {-^ [y - y*{t)] M;' [y - y.(t)] 



where y*(y', t) is the mean value of the phase space vector y and Ai^{t) is the variance matrix. 
The explicit form of y,,, and Ai^ can be obtained by solving the equations (which are immediate 
consequences of the Fokker-Planck equation ( |A.1|) ) 

M, = AM, + MJl + 2B (A.4) 
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y = Ay + K (A.5) 

with the initial conditions A^^,(0) = and y*(y', 0) = y', which are equivalent to the property 
G (y, y', 0) = 5 (y — y') distinguishing the propagator from all other Gaussians. 

In the case under study equations (|A.4|) coincide with the set (|42|) -(^5 D . Putting their 
solutions (^BD-(|5DD to the right-hand side of ( |A.^j| ) we obtain the propagator (|TBp. 
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Figure 1: The 'purity' jj, and the 'coherence' C versus the 'compact time' u=l — e"^'*'* for the 
initial coherent state with jap = 1 and for three diferent temperatures. The order of the curves 
from top to bottom (in the left-hand side of the figure) : C for u ~ (zero temperature) (in this 
case 11= 1); II for = 1, C for i/ = 1; for i/ = 10, C for u — 10. 



Figure 2: The 'purity' and the 'coherence' C versus the 'compact time' u = 1 — e~ for 
the initial coherent state with lap = 10 and for three diferent temperatures. The order of the 
curves from top to bottom (in the left-hand side of the figure): C ior v — (zero temperature) 
(in this case = 1); for i/ = 1, C for = 1; /x for i/ = 10, C ior u — 10. 



Figure 3: The 'purity' /i and the 'coherence' C versus the 'compact time' u = 1 — e"^'''* for the 
initial odd coherent state {(f — tt) for different values of the parameters a = lap and i/. The 
order of the curves from top to bottom (in the right-hand side of the figure) is as follows. I - 
// for a = 1 and z/ = 0; II - /i for a = 10 and u = 0; III - C for a = 10 and u — 0; IV - C for 
a — 1 and u = 0; Y - /j, ior a — 2 and = 5; VI - C for a = 2 and u — 5. 



Figure 4: The 'accompanying coherence' ^ versus the 'compact time' u = 1 — e for the 
initial even coherent state (</? = 0) with \a\'^ — 10, for three different temperatures: u — 0, 1, 10. 



Figure 5: The 'purity' n and the 'coherence' C versus the 'compact time' u = 1 — e^'^"'^ for 
the initial squeezed coherent state with p — \a\'^ — 1 and — 7r/2, for two values of the 
equilibrium mean photon number in the reservoir: z/ = (two close curves relatively far off the 
bottom) and z/ = 2 (two close curves nearby the bottom). In each pair of close curves the upper 
one corresponds to /i while the lower one gives C. The curves inside the internal box give the 
same functions versus the usual (scaled) time r = 27^. 



Figure 6: The 'thermalization coefficient' V versus the 'compact time' u = 1 — e"^'"'* for the 
initial odd coherent state ((/? = vr) with different values of the parameter \a\'^ — 1, 20; in the 
low temperature {y = 0.01) and high temperature = 10) cases. 



Figure 7: The 'thermalization coefficient' V versus the 'compact time' u = 1 — e"^'"'* for the 
initial squeezed coherent state with p = 3, lap = 1 and = 7r/2, for two values of the 
equilibrium mean photon number in the reservoir: v — 0.01 (lower curves) and v = 2 (upper 
curves). The curves inside the internal box give the same functions versus the usual (scaled) 
time r = 27^. 
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Figure 8: The 'thermalization coefficient' V versus the 'compact time' u = 1 — e~^'^* for the 
initial coherent and Fock states with equal mean numbers of photons: M = |ap = 1 and 
M = lap = 20, in the low temperature {v — 0.01) and high temperature {v — 10) cases. I - 
the Fock state with M = 1 and v = 10, II - the Fock state with M = 20 and v = 10, III - the 
Fock state with M = 20 and v = 0.01, IV - the coherent state with |ap = 20 and z/ = 10, V - 
the Fock state with M = 1 and u = 0.01, VI - the coherent state with |ap = 1 and u = 0.01, 
VII - the coherent state with |ap = 20 and v = 0.01. The internal box shows the behaviour of 
the last two curves in another scale. 
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